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1 Introduction 



Recently Dijkgraaf-Vafa proposed that a holomorphic information in M = 1 gauge theories 
with classical gauge groups derive from matrix model [H El El Hj • In particular the effective 
superpotentials for the gauge theories are given by the genus zero free energies which are 
the summations of planar diagrams in matrix models. This perturbative sum of planar 
diagrams turns out to be a nonperturbative sum over fractional instantons in the M = 1 
gauge theory. 

Dijkgraaf and Vafa have reached this duality via string theory route using the model 
with one adjoint matter. It was discussed that in [3J E] large N dual of U (N) Af = 2 
gauge theory deformed by certain tree level superpotential is realized as type IIB string 
theory on Calabi-Yau threefold with fluxes. The discussion for type IIB superstring was 
extended in [HI El EE EE HH CEZ1 UHj- The equivalence of effective superpotential for 
Af = 1 gauge theory with that of type IIB superstring theory on Calabi-Yau manifold 
with fluxes [13 d UHJ was proved in [TTHTH]. 

There are many discussions on this D-V conjecture [19]- [73] . The generalization of 
D-V duality to the gauge theories with massive flavor was discussed in jHHl EH EH EH HH 
ESI EH ESI EZl 1121 ESI EH! EH1 - In JSHIEOI perturbative analysis of the matrix model, which 
corresponds to U(N) gauge theories with arbitrary tree level superpotentials. Af = 2 
information was derived these perturbative results including the Seiberg-Witten curve. 
The analysis beyond the planar limit was discussed in [221 EH]- in this note we apply 
these perturbative result with SO/Sp gauge theories. 

In SO/Sp gauge theories were studied using the superspace method that was de- 
veloped in [33] . In [331 EH] antisymmetric matrix model was considered as a generalization 
of DV to SO/Sp gauge theory. On the other hand real symmetric matrix model was used 

in [22EEI- 

The organization of this paper is as follows: In section 2, we review the discussion 
in [HSj and generalize to the case with massive flavor. In section 3 we give perturbative 
result for the effective superpotential of SO/Sp gauge theory in terms of the relations 
discussed in section 2. In section 4 we give exact solution for the case with quartic tree 
level superpotential with massive matter. In section 5 from the matrix context we derive 
Seiberg-Witten curves for the gauge theories. 

Note added: After completion of this note, we have received [20], which have some overlaps 
this paper. 
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2 Relations among Free Energies 



In this section we represent some relations between free energies. These relations are used 
following section for the computation of effective superpotentials. 



2.1 General Property 

In the standard large N expansion the partition function of matrix model has topological 
expansion, 

Z = exp (5>7*F X ) , (2.1) 

where x is Euler number of the surface. In this letter we will take into account the surface 
with crosscap, because we discuss SO/Sp gauge theory. Euler number with crosscap c is 
described as x — 2 — 2g — b — c, where g, b, c are the number of genus, boundaries and 
crosscaps respectively. The leading contribution comes from S 2 with Euler number x — 2. 
We denote this free energy as JF . The Euler number of RP 2 is x = 1 that is equal to 
the Euler number of orientable one boundary surface, disk topology. We denote this one 
crosscap contribution and one boundary contribution as 

F X=1 = Q + J 7 1 , Go = 9s1ogZ\ RP 2, T x = g s \ogZ\ disk . (2.2) 

Now let us define matrix models corresponding to SO/Sp gauge theories. The partition 
function of matrix models is given by, 



1 1 Nf 

= /rf$exp Wo(^)--J2(Qi^Q I + m I Q I Q I ) 

J 9s 9s j = i 



(2.3) 



where $ is 2M x 2M real antisymmetric matrix 1 or real matrix with relation, $ r 
(J$J) nm where J denotes the Sp invariant skew symmetric form 



Imxm 
— Imxm 



J 2 — ~^2Mx2M- (2.4) 



Q, Q is 2M dimensional vector. For simplicity we dropped volume factor for the gauge 
group, which gives Veneziano-Yankilowicz term. We assume that tree level superpotential 
Wo is even function, i.e. 

n+l 

W o = E f£ Tr$2fe - (2-5) 



1 Note that matrix size 2M have not relation to gauge group rank. 
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On these two matrix models and Hermitean matrix model that correspond to U(N) 
gauge theory, some relations are already known. For the free energy coming from disk 
topology 

F 1 so (g)=^ p (g) = ^(g' = 2g). (2.6) 

where gauge groups denote the matrix models corrsponding to these gauge theories. The 
argument g,g' denote coefficients of tree level superpotential. Note that for SO/Sp case 
the tree level superpotentials are even function then using this relation we assume that 
the matrix model corresponding to U(N) gauge theory also have even function tree level 
superpotential. And for genus zero free energy we have 

^°(f) = Fo P U) = \^{9 = 2/)- (2.7) 

Multipicating 1/2 over all factor to and replaceing Wq to 2Wq, namely multipling 2 
to all coefficent of superpotential, it equals to Tq° and J 7 p '. Next relation is the one 
between JF and Go, 

where — sign for antisymmetric matrix and + sign for the matrix corresponding to Sp 
gauge theory respectively and So is 't Ho oft coupling that defined next subsection. These 
relations were checked in j33] perturbatively in the viewpoint of field theory using su- 
perspace formalism and in (23 EH] using loop equation for antisymmetric matrix models. 
We give another derivation for these relations. In [SB] the relation between the T§ corre- 
sponding to the U(N) gauge theory and the one SO(N) gauge theory. 

With these relations we can obtain effective superpotentials for SO/Sp gauge theory 
from matix models. In EH] the experssion was given, 

aT so/s P r, T so/s P FlT so/Sp 

W« = E N^h— + 4g ° + K = * ^^k— + E N ^hr~ + (2-9) 

i=0 "^0 i=i CDi 

It is interesting that two free energy Q and T\ that are the same Euler number, appear 
to effective superpotential having different coefficient. For the latter convenience let us 
introduce new notation iVj, which means that N = N =p 2 and iVj = iVj i = 1, • • •. In 
section 3 and section 4 we will see effective superpotentials for multicut case and one cut 
case, using these formulae. 



3 



2.2 Matrix Model for SO(N) with Flavor 



In this subsection we discuss the antisymmetric matrix model. These matrix model was 
studied in jHHl IHEj f° r the generalization of DV proposal to the SO/Sp gauge theory. 
We want to extend these discussion to the model with massive matter. The partition 
functions of these matrix models are given in ()2.3|) . As discussed in [33] we can integrate 
out Q and diagonalize matrix $, 



$ -> UDU 1 



-D 



2k - 1, k 



N) others zero.(2.10) 



where we used notation ±A,-, % 



N as eigenvalues of matrix $. Then A, are pure 



imagenary value. Using these variables Aj we can rewrite as follows, 



M 



iy ($)^E 2 ^o(A fc ), d$ 
k=i 



M 
k<l 



M 

n<*Ai. 



(2.11) 



. M ( 9 M N f M \ 

Z= /n^exp --E^(A fc )+2Elog(A 2 fe -A 2 )-EE lo g(A 2 -m 2 ) . (2.12) 

J i \ »* fc=l fc<Z 7=1 i / 

where third term come from the factor det _1 ($ + rrij). The equation of motion for Aj is 
given by 

- ^o(Aj) + 2 E t^t, " E tF^ = 0. (2.13) 

fl's jjii \ - A j i=i A i ~ m I 

As usual let us introduce resolvent u, 



1 _ / 1 \ x M 



UJ = 



Tr - ^ = T7 ^ \2' (2^) 



2M " Vx-$7 M^x 2 -A, 2 



In the large M limit with g s M = S fixed, it is convenient to introduce density of eigen- 
values, p(A) = ]gSj5(iA — i\i),J p(X)d\ = 1, where i is necesary in the delta function 
becauce Aj are pure imaginary. With this new variables we can write loop equation as 
follow, 

_ -^(x)^(x) + — /(x) = 0, f(x) = — E l A 2 + x2 1 % (2-15) 

where third term in ()2.13j) supressed in this limit. As in [JJ], introducing new variable, 
y = Su(x) + Wq(x), we can rewrite loop equation as hyperelliptic curve, 

y 2 = Wg l -f(x). (2.16) 
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Generally this Riemann surface has 2n + 1 cuts, we can express the filling number of 
eigenvalues in terms of the integral around the cut, 

M = m/ ydx, Si = g B Mi, S = ^M . (2.17) 

Ji th cut Z 

Since Mi's have 2M = M + 2£;M;, these Si 8 have the relation, S = S + J2iSi- We 
can rewrite partition function up to constant term, 



Z = exp ( J d\p(X)W (X) + ^ J dXdpp(X)p(p) log(A 2 - p 2 

N i 

Y J dXp(X) log(A 2 - mj) (2.18) 



9 S M N ' 



9s I=1 

Then we can read off the contribution from disk from ()2.18|) as follows 

S 

i=i 

We can rewrite this integral as follow 

Nf r rP 2r 9 Nf 

dXp(X) 

J Jmj A — X Z 



a N f 

K = "f E J ^V(A) log(A 2 - mj). (2.19) 



*i = -|E/^P(A)/ 3^-,-|E/^p(A)log(A 2 -P 2 ) 

= / dxu(x) - -N f D(P) (2.20) 

J— mi Z 

= -E / Z/^ + -^iy (P)-E^o(m / )--iV /J D(P) (2.21) 

where we defined function D(P) = J dXp(X) log(A 2 — P 2 ) and P is one of two infinity on 
the x-plane. 



3 Perturbative Analysis 

In this section we give effective superpotential for the SO/Sp gauge theories with arbitrary 
tree level superpotential up to first instanton correction. In the previous section we study 
the relation between real symmetric matrix model and unitary matrix model. Using the 
relation we can obtain real symmetric matrix model result form unitary matrix model 
result and explicitly give effective superpotential for the SO/Sp gauge theories. Since in 
[571] perturbative expression for the genus zero free energy and disk amplitude of unitary 
matrix model was given, we can give the ones for real symmetric matrix model. For this 



5 



purpose we rewrite equations in [SU] for the case with even tree level superpotential and 
symmetric distribution of eigenvalues. We write the tree level superpotential as 

n+l 



92k rp„^2fc 

2k 

Let us introduce new variables e 



Wb(*) = E ^Tr$ 2 *, W&x) = & +2 xE(x 2 - al). 
k=i ZK k=\ 



(3.1) 



e = 0, ti = di 



-Q{, i = 1, 2, • • • n. 



(3.2) 



We evaluate matrix integral perturbatively about and extremal point $ = $o ; Qo 
0,Q = 0. 

/e_„l Mn ••• \ 

e_„ + ilM n _i ••• 



V o 







enlM n / 



S* = g- s Mj = S- h g s 0, 5j = fixed. 



(3.3) 



(3.4) 



We fix the gauge ^ = (i / j)and introduce Grassmann-odd ghost matrices B^ and 
Cij. The partition function is given by the gauge-fixed integral 



Z = * exp ( -— W ($o) ) / d^udBijdCijdQ^QjexpiSuin + S int ) 
volyLr) \ q s I J 



(3.5) 



kin 



#2n+2 



n Nj 



E ^ Tr ^ -EE %Tr(%Cy --EE fuQuQi » (3-6) 



mt 



n 2n+2 



-n /=1 



i=— n p=3 



V 



n I n Nf 

~ E Tr^^Cy - BjiQj^jj) - - E 

i=— n 9s i = — n j 



where Ri = Yij^i e ij with = e« — e^, /*/ = e« + mj and 



7p,i 



(p-1)! 



'A 
dx 



p-1 



n ( x - e fc) 



k=—n 



(3.7) 



(3.f 



Note that these relations are the one for the unitary matrix model. Then the effective 
superpotential given by this matrix model is the one for U(N) gauge theory. However 
we can obtain real symmetric matrix model free energy replacing #2n+2 to 2g2n+2 and 
multiplying the over all factor 1/2. 
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3.1 Genus zero free energy 



Using the relation between Fq U and discussed in section 2, we can get the following 
result, 



(3) 



1 



1 1 n ^ 3 / i \ i n c3 



i=—n kj^i 



q2 C 1 n C C C 

J i "Ji; ^ _|_ 2 5^5^ °i°kOl 

Ri&ik i^-i &il i=—n. hJzi lJii Ri&ik&il 



i=—n k=£i l^i 



EE 

i=—n k^i 



C2 C 
Ok 

Ri e ik 



(3.9) 



where we denoted free energy of third order in 5* as JFq 3 -* . These result comes from the 
following Feynman diagrams. 



+ 






+ 




+ 




+ 




Figure 1: 2-loop diagrams contributing to the third order of S. Striped double lines correspond to ghost 
propagators. Note that we does not compute unoriented diagram but modify the result for Hermitean matrix 
model with even function tree level superpotential. 



Though we are considering the tree level superpotential, which is even functions poly- 
nomial, perturbative expansion around classical solution gives potential, which is not even 
function, (|3.7|) . Thus we have cubic interactions in Feynman diagrams. 

We can check this result for the special case, i.e. <3> 4 case. Let 0, ±ai are the critical 
point for the tree level superpotential. Thus e_i = — a 1; e_ = 0, e x = a x and = 
2ctx, Rq = — a 2 , R\ = 2a 2 , 



•r 



(3) 



+ 




af 



C C2 C-i 

af 2af 



~Sq — 2,SqS^ -\- 45'g>S'i 




4g,af V 2~ u ~ u " 1 ' ~ "7 " (3 ' 10) 
We will see that this result agree with the one in [5] obtained by Calabi-Yau with fluxes. 
We review the result in Appendix of computation of period of Calabi-Yau manifold 
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with fluxes. The computation was done in semi classical limit A<1. Then the period of 
compact cycle that interpreted as glueball superfield was small. The period of noncompact 
cycle was described as the perturbative expansion of the period of compact cycle up to 
4-th order. Using the prepotential perturbative part of the result is written as follows, 



OS, 



(3.11) 



V 



CY 



4gaf 



2 



2S S( + 45^1 + 



8g 2 af V8 



S% - US^S 1 + 185 2 5i - ASoSf 



+ 



16# 3 a} 2 



9 c5 , ^33 4 304 2 3 524 3 2 40 4 , 

"2 ^"^O ! "1 ^"^O ! ""^"^O ! "^" D 0'Jl J 



(3.12) 



First order terms agree with (|3.10|) . In [5] using these perturbative expansion of periods 
effective superpotential coming from the flux was given. Then integrating out glueball 
superfield effective superpotential that was written as A and parameters in tree level 
superpotential was described. These superpotential can be compared with confining phase 
superpotential given purely field theory analysis. A few ensamples show agreement these 
two result. Then the analysis of matrix model describes precisely field theory result. 

3.2 Disk amplitude 



For the disk amplitude there is no difference for two models, then the perturbative result 
is given in terms of the results in [50J, 



TV 



(2) 



N 



f 



2 EE^ 

1=1 jjki ih i c -ijjti 



1=1 \i=- n -L^JH j^i ^13 

These result comes from the following Feynman diagrams, 



1 n q2 

2 i= _ n Rifip 






(3.13) 



Figure 2: 2-loop diagrams with one boundary contributing to the second order of S. Fat lines represent a 
boundary that comes from the propagation of Q's. As in the case we consider Hermite matrix model and 
then modify the result for the real symmetric matrix model case. 



Let us also for this amplitude give explicit representation for the $ 4 case, 



(2) 



Nf r 

E 

1=1 L 



C<2 
D 



2a 2 m 2 



+ 



2SqSi 



a 2 (m 2 — a 2 ] 
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+ 



S 2 



m 



2\2 



(3.14) 



Finally we obtain effective superpotential from matrix model calculation up to second 
order in Si, namely up to first instanton correction, 

n 

w eS = J2Si 

i 

where first term comes from volume factor |67j . 




+ 



dS,, 



(3) 



+ -F 1 (2) + -- 



(3.15) 



4 Exact Effective Superpotential for Unbroken Gauge 
Group 

In this section we discuss one-cut solution for SO/Sp case with flavor. For $ 2 case there 
are exact one-cut solutions in O HZl EH1 • We want to generalize these results to a quartic 
tree potential case, 

W = — $ 2 + -$ 4 . (4.1) 
2 4 y J 

The one-cut Riemann surface is already discussed in P, 



y = y/w? + f{x) = g{x 2 + A + 2 / x 2 ) v /(x - 2n)(x + 2/i), (4.2) 

where A = One of two cut shrinks to be singularity. Using this y we can obtain 
following relations. It behaves as 1/A when |A| goes to infinity. This constraint gives 
following relation, 



S = 6gf + 2M/i 2 , ^ = " f + 1 + ^ - jHS)- (4-3) 
With the relation given in section 121 we can calculate the effective superpotential 



dS hp. 



/ dx(x 2 + A + 2/i 2 ), x 2 - 4/i 2 - W (P) + -D(P) 

J2u 2 



-S log A + | log ^ - j - 9 -^-L(S) + | log ^L(S) + ^D(P), (4.4) 



Nf Mo N f 1 N N 9 

Pi = J2 dx(x 2 + A + 2/i 2 ) ^ - 4// 2 - £ r^oK) + -ttW(P) - -^-D{P) 



X; r m/ 2 2 3 3/i 2 + A 



E 

i=i 



I v m 2 — 4/i 2 ) 2 + — ' rrij^m 2 — 4/i 2 
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-2/i 2 (3/i 2 + A) (bg2A - log(m 7 - 



m 



2 - V 



-W(mi) - -D(P) 



N f r 

E 

1=1 L 



m 
4 



(4.5) 



A 2 / 
- T L(S)(l 



L(sy 



log 



2A 



l -W(mj) - ^D(P) 



where if (5) = Jl - ^L(S) 



Since logA divergent piece can be renormalized to bare coupling a, we can replace A 
to A, which is interpreted as dynamically generated energy scale. Until now we used A 

3 

for simplicity, but correct dimesional parameter is A 2 so we replace the parameter. And 
using the matching relation, A 3N = det mA 3N ~ N f , where m is massmatrix that eigenvalues 
are mj, we can obtain effective superpotential, 



W, 



eft- 



s' log 



A^ 

S n 



+ N 



T -^(S) + -log-L(5) 



Nf 



+ E 



^K(Sf + ^ (l + A- (S) 

^ ) ( l0g (T^(S))) 



A 2 / 

- T ^(S) 1 



(4.6) 
(4.7) 



Now let us consider g — > limit. There are many discussions for this model with U(N) 
gauge theory. In this limit we obtain the following result. The effective superpotential 
for this case is described as 



W. 



eft 



Slog 



A 



3N 



S N 



+ NS 



(4.t 



N f 



^ s 



Mm] 



'1- 



45 



Mm? m /l 1 r 4S 



4 V Mm? 

This is exactly the same as the result from the discussion for flux with Calabi-Yau ge- 
ometry in [IB]. And If we have introduce new function 2fi(S) = 1 
rewrite as, 



4S 
Mm! 



we can 



W, 



eft 



Slogf^l +NS-Y (--—?— -logfAS)) 

I s J 1 V 2 2 //(s) 7 



This result agree with matrix model analysis in 



(4.9) 

up to non-ecessial convension of A. 
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5 SW Curves from Matrix Models 



In section 2 we reviewed the knowledge that for introducing the contribution RP 2 we 
have only to replace iV to N =p 2. Thus with this lesson we can straight forwardly extend 
the derivation for Seiberg-Witten curve for the U(N) in [5U| to SO/Sp model. To get 
M = 2 information, at first we consider M = 1 theory with tree level superpotential is 
W (x) = Y,i N+2 Tr$ 2fe and then choose the special vacuume, N { = 1, i — 0, 1, • • • , N. 
Next we should turn off the tree level superpotential, namely g2N+2 — > 0. In this vacuum 
effective superpotential is given by 



N Q-TT N 



W eS = (l T 2)|g+E^-27r^ES 

N l N 

= (1=f2)/ ydx + J2 ydx--ToJ2 y dx + const, (5.1) 

where Bi is the closed loop through z-th and iV-th cut. 

As in U(N) case we use equation of motion for &2i's % = 1, • • • , N. The derivative of y 
with respect to £» 2 i is holomorphic function and form a complete basis. Let us introduce 
one forms, which form complete basis to holomorphic differentials, 

° !i dx = -—x 2n dx = Cn- (5.2) 



db 2n 2y 

Using this notation the equation of motions, dW^jdbm = 0, can be written as 

„ rQ A rP Nf r-mi 

N ( k -(N-N f ) ( k - / Cfc = 0, (mod period integral). (5.3) 



Thus from Abel's theorem there exists a function z(x) on the Riemann surface with an 
th order pole at Q, an (N — Nf) th zero (or pole if Nf > N) at P, and simple zeros at 
—mi. For the N > 2Nf we find the function z(x) as follow, 

z(x) = A{x) - ^(A(x)) 2 - Cdet 2Nf (x + m), for SO(2N) (5.4) 
z(x) = B(x) - ^B(x) 2 - Ddet 2Nf (x + m), for Sp(2N) (5.5) 

where A(x), B(x) are airbitrary polynomials defined below and C, D are constant. These 
function have N th order pole at Q and — 2Nf th order zero at P and 1 th order zero 
at —mi respectively. For z to be a function on the Riemann surface the square root in z 
must be proportional to y, 



y 2 
y 2 



A(x) 2 - Cdet 2Nf (x + m)\ T , (5.6) 
B(x) 2 - Ddet 2Nf (x + to)] (T Sp ) _1 , (5.7) 
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where T's are the proportional function. 

We can determine these quantity by taking classical limit A —>■ 0. In this result matrix 
model curve y becomes Wq. So It is natural that A(x) and B(x) are propotional to 
charactristic function P 2 n, which becomes Wq under the classical limit, 

Afar) = A?, B(x) = (x 2 P 2N (x) + 2A m+4 ~ 2N fPim). (5.8) 

x 2 

For the constants C, D classical value of Qi, Qi are zero, it is natural that C, D are 
proportional to A. With these relation we have the following results in this limit, 

T SO = T S P = x 2_ (5 g) 

Then we have the following result, 

y 2 = P 2N (xf - x 4 A bo det 2Nf (x + m), (5.10) 
x 2 y 2 = (x 2 P 2N (x) + 2A 26o Pfm) 2 - A 26o det 27V/ (a; + m). (5.11) 

These are the Seiberg-Witten curve for the gauge theory. Here we have derived these 
results using only matrix-model methods. 
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